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^ ■ Abstract 

We present a variational basis-set calculational scheme for elastic scat- 

(N . 

^ . tering of positronium atom by helium atom in S wave and apply it to the 

. calculation of the scattering length. Highly correlated trial functions with ap- 

. propriate symmetry are used in this calculation. We report numerical result 

for the scattering length in atomic unit: (1.0 it 0.1)ao. This corresponds to a 

O , zero-energy elastic cross section of (4.0 it 0.8)7raQ. 
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I. INTRODUCTION 



Recent successful measurements of ortho positronium (Ps) scattering cross sections by 
H2, N2, He, Ne, Ar, C4H10, and C5H12 [|l|-0| have spurred renewed theoretical activity in this 
subject iP Hl^j . Of these, the Ps-He system is of special interest as it is the simplest system 
in which there are experimental results for total cross section and pickoff quenching 
rate . The experimental results for partial and differential cross sections for this system 
should be available soon A complete understanding of this system is necessary before a 
venture to more complex targets. 

The pioneering calculations in this system using the static exchange approximation were 
performed by Barker and Bransden p^JT6|l and by Fraser and Kraidy |l7,18]. There have 



also been R-matrix |Q, close-coupling (CC) |1TT|JT^ and model-potential |19] calculations for 
Ps-He scattering. More recently, there has been successful calculation of Ps scattering by H 
0, He [gTHH' Ne Ar and H2 using a regularized model exchange potential 
in a coupled-channel formulation. 

However, there is considerable discrepancy among the different theoretical Ps-He cross 
sections at zero energy which we discuss below. The static-exchange calculation by Sarkar 
and Ghosh and by Blackwood et al. g] yielded 14.387rao (at 0.068 eV), and 14.587rao 
(at eV), respectively, for the elastic cross section. The inclusion of more states of Ps in 
the CC and R-matrix calculations does not change these results substantially. The 
pioneering static-exchange calculations by Barker and Bransden yielded 13.04 ttoq and 
by Fraser ||16| yielded 14.27raQ for zero-energy Ps-He cross section. These results are in good 



''0 

agreement with each other. However, the model potential calculation by Drachman and 
Houston [1^ yielded T.TSttoq and by this author yielded 3.347raQ for the zero-energy Ps- 
He cross section. So there is considerable discrepancy in the results of different theoretical 
calculation of low-energy Ps-He elastic scattering. 

On the experimental front, there have been conflicting results for the low-energy Ps-He 



elastic cross section by Nagashima et al. [Q, who measured a cross section of (13 ± 4) 



at 0.15 eV, by Coleman et al. j^, who reported 97raQ at eV, by Canter et al. who 
found 8.477raQ at eV, and by Skalsey et al. who measured (2.6 ± 0.5)7raQ at 0.9 eV. It 
is unhkely that these findings could be consistent with each other. 

The results for the total cross section of Ps scattering obtained from the coupled-channel 
calculation employing the model potential ||2l|-p3| are in agreement with experiments of Refs. 

at low energies. For Ps-He, this model, while agrees pTl-[23[] with the experimental 
total cross sections in the energy range to 70 eV, reproduces ||2^ successfully the 
experimental pickoff quenching rate |]5|,p!^. All other calculations could not reproduce the 
general trend of cross sections of Ps-He scattering in the energy range to 70 eV and 
yielded a much too small quenching rate at thermal energies |1^,|1^J2^. However, the very 



low-energy elastic cross sections of the model-potential calculation pl|-p^ are at variance 
with the experiments of Refs. . 

Pointing at the discrepancy above among different theoretical and experimental studies, 
Blackwood et al. called for a "fully fledged calculation" to resolve the situation. Here 
we present a variational basis-set calculational scheme for low-energy Ps-He scattering in 
S wave below the lowest Ps-excitation threshold at 5.1 eV. Using this method we report 
numerical result for the scattering length of Ps-He using a one parameter uncorrelated He 



ground-state wave function p6 |. 

We present the formulation for the variational basis-set calculation in Sec. H, the nu- 
merical result for Ps-He scattering length in Sec. HI and a summary in Sec. IV. 



II. FORMULATION 

Because of the existence of three identical fermions (electrons) in the Ps-He system, one 
needs to antisymmetrize the full wave function. The position vectors of the electrons — ri 
of Ps, and and of He — and positron (x) measured with respect to (w.r.t.) the massive 
alpha particle at the origin are shown in Fig. 1. In this configuration the wave function for 
elastic scattering in the electronic doublet state of Ps-He is taken as 



^^(1, 2, 3) = 0(2, 3)r7(l)Fk(l, 2, 3)x(l, 2, 3) (1) 

where k is the incident Ps momentum and 

2, 3) = a(l) X -L[«(2)/5(3) - /5(2)a(3)], (2) 

represents the doublet wave function of Ps-He and where a denotes spin up state and [3 down 
and 7/(1) denotes the Ps wave function of electron 1. The He ground state wave function 
0(2, 3) and the scattering function -Fk(l, 2, 3) are symmetric under the exchange of electrons 
2 and 3. The spin function x(l, 2, 3) is antisymmetric under the same exchange. The full 
antisymmetrization operator for the three electrons is (1 — P12 — -Pi3)(l — -P23) where Pij 
is an operator for exchange in both space and spin of electrons i and j. As the scattering 
wave function (|I]) above is already antisymmetrized with respect to electrons 2 and 3, the 
operator (1 — P23) in the antisymmetrizer is redundant and the relevant antisymmetrizer in 
this case is ^ = (1 — P12 — Pis). Hence, the fully antisymmetric state ip^ of Ps-He scattering 
is given by 

^if = Ml = 0(2, 3)r7(l)Pk(l, 2, 3)x(l, 2, 3) - 0(1, ?>)v{2)F^{2, 1, 3)x(2, 1, 3) 

-0(l,2)r/(3)Pk(3,2,l)x(3,2,l). (3) 

The projection of the Schrodinger equation on the doublet state ^(l, 2, 3) is 

(X(l, 2, 3)|(i/ - E)\i^:^) = {H- E)|0(2, 3)r7(l)Pk(l, 2, 3) - 0(1, 3)r/(2)Pk(2, 1, 3)) = 0, (4) 

with H the full Ps-He Hamiltonian. The incident Ps energy E = 6.8A;^ eV. Using the 
identities (x(l, 2, 3)|x(l, 2, 3)) = 1 and (x(l, 2, 3)|x(2, 1, 3)) = (x(l, 2, 3)|x(3, 2, 1)) = 1/2, 
we see that the two terms on the right-hand side of Eq. (^ give equivalent contribution 
which are combined in Eq. (^ which is rewritten as 

-E){1- P12) 10(2, 3)r/(l)Pk(l, 2, 3)) = 0. (5) 

Hence after the spin projection to the doublet state the effective antisymmetrizer to be used 
on state (|l]) is Ai = (1 — P12). We shall use this antisymmetrizer in the following and supress 
the spin functions. 



The full Ps-He Hamiltonian H can be broken in the convenient form as follows: H = 
Hi + Vi where Hi includes the full kinetic energy and intracluster interaction of He and 
Ps for the arrangement shown in Fig. 1 and Vi is the sum of the intercluster interaction 
between He and Ps in the same configuration: 



Vi 



2 2 1111 

X n ri2 P2 ri3 Ps 



(6) 



We employ the position vectors Sj 



(x + rj)/2, pj = x-r^, = - r^, i,j = 1,2,3, 



The fully antisymmetric state satisfies the Lippmann-Schwinger equation |27 



K) = m+GiMi\iJi), 



(7) 



Ml = ViAi + {E~ Hi){l - Ai) = AiVi + (1 - Ai){E - H 



1), 



where the channel Green's function is given by Gi = {E + iO — Hi)^^ and the incident wave 
satisfies {E — Hi)\(f)l^) = 0. We are using atomic units (au) in which ao = e = m = h = l, 
where e (m) is the electronic charge (mass) and oq the Bohr radius. 

The properly symmetrized transition matrix for elastic scattering is defined by 
{K\T^\K) = (KlVil^k) = (0kl^iA|O = (V'klA'^il^) H- A basis-set calculational 



scheme for the transition matrix can be obtained from the following expression p8 



''^T^m = {iji\AiVi\<Pi) + {<Pl\AiVi\i^l) - ki^i\AiVi - MiGiAiVMi). 



(9) 



Using Eq. (|^), it can be verified that Eq. (j^) is an identity if exact scattering wave functions 
are used. If approximate wave functions are used, expression (j^) is stationary w.r.t. small 
variations of l^/^^) but not of This one-sided variational property emerges because of the 
lack of symmetry of the formulation in the presence of explicit antisymmetrization operator 
A\. However, this variational property can be used to formulate a basis-set calculational 
scheme with the following trial functions 



N N 



(10) 



where the suffix t denotes trial and fn,n = 1,2, ...,N, are the basis functions. Substituting 
Eq. (|10D into Eq. (|) and using this variational property w.r.t. we obtain ||2^ 



N 



m=l 



{D-')mn = {fmlAiVi - [A^i + (1 " Ai){E - H,)]G^AlVi\fn). (12) 

Using the variational form ([Til) and definition {4>k\T-^\(j)l^) = (^/'kl-^i^il^k) obtain the 
following basis-set calculational scheme for the transition matrix 

N 

''^T-^\K)t= E {K\AiV,\U)D^^{UA,V,\^l). (13) 

m,n=l 



Eqs. ( |12|) and ([13|) are also valid in partial-wave form. 



In the present S-wave calculation, the basis functions are taken in the following form 

/m(r2, rs, Pi, si) = (p{r3)gm{r2, pi, Si), (14) 
^m(r2,pi,si) = y?(r2)r7(pi)e~^-"^-"-^^~^-^^~^-(^^+^'^^)~^-("+'^^) sm(fcgi) ^ 

where 5^, ^m, /^m, 7m, and /i^ are nonlinear variational parameters. The ground-state wave 
function of the He atom is taken to be 0(r2, r^) = (/?(r2)(y9(r3) with (f{r) = A^^^ exp(— Ar)/-\/7r 



and A = 1.6875 and r](p) = exp(— 0.5p)/v Stt represent the Ps(ls) wave function. For 
elastic scattering the direct Born amplitude is zero and the exchange correlation dominates 
scattering. To be consistent with this, the direct terms in the form factors (/m|-^i^i|0k) 
and (0k|^iVi|/„) are zero with the above choice of correlations in the basis functions via 
7m and fim- This property follows as the above function is invariant w.r.t. the interchange 
of X and ri whereas the remaining part of the integrand in the direct terms changes sign 
under this transformation. In Ps-He elastic scattering the electron 2 of He is the active 
electron undergoing exchange with the electron 1 of Ps whereas the electron 3 of He is a 



passive spectator. In this calculation we include in Eq. (|T5D correlation between electrons 
1 and 2. Consequently, we deal with integrals in three vector variables - ri, r2 and x. If we 
also include correlation involving electron 3 we shall have to deal with integration in four 



vector variables, which is beyond the scope of the present study. However, we beheve that a 
meaningful calculation can be performed only with correlation between the active electrons 
1 and 2. Hence, to avoid complication we ignore correlation involving electron 3, which is 
expected to lead to correction over the present study. 

In S wave at zero energy, sin(A;si)/(A;si) = 1 in Eq. (|T5|); also, = 
ip{r2)(p{r3)ri{pi)sm{psi) / (psi) . The useful matrix elements of the present approach are ex- 
plicitly written as p8[ 



1 f S1T\.(t)S2) 

{(l)l\AiVi\fn) = -— / '^(ri)<^(r3)r/(p2) [\^i]/„(r2, rg, pi, s^drac/rgrfpirfsi, (16) 

^ ZTT J pS2 

= I '/5(ri)r7(p2)^^^^^^^[Vi]5f„(r2,pi,Si)dr2dpirfsi, (17) 

ZTT J pS2 

{U\A,V,\<f)l) = I (7^(ri,p2,S2)[Vi]y.(r2)r7(pi)^^^^«pif/si, (18) 

^ In J psi 

{fm\AiVi\fn) = J gmiri, p2,S2)[Vi]gn{r2, pi,si)dr2dpidsi, (19) 



with 

Vi 



h{x) - h{ri) + — - — 
ri2 P2 



h{x) = -+ "^""P^ ^^"^^ + A exp(-2Aa;), (20) 

X X 



2 

ifmmCAiV.lU) ^ — dp{U\AiVi\<Pl){<Pl\AiVi\fn), (21) 

where the so called off-shell term (1 — Ai){E — Hi) has been neglected for numerical sim- 
plification in this calculation. This term is expected to contribute to refinement over the 
present calculation. In this convention the on-shell t-matrix element at zero energy is the 
scattering length: a = (0g|T-^|0g). 



All the matrix elements above can be evaluated by a method presented in Ref. p9 
We describe it in the following for {(f)p\AiVi\fn) of Eq. (|1^). By a transformation of vari- 
ables from (r2,pi,Si) to (si,S2,x) with Jacobian 2^ and separating the radial and angular 
integrations, the form factor (|T6|) is given by 

26_x3 „oo ^sm{ps 



^Jl^i^il/n) = rds2sf-^^ Hdsisle-^-^^ Hdxx'e-^-^ 
^ lon-^ Jo ps2 Jo Jo 

X y"e-('"^^+'"'^/2)g-(cr2+dp2/2)g-7„ri2[y^|^^^^^^^^^ ^22) 



where a = A + /i„, b = 2a„ + 1, c = A + 5„ and d = 27„ + 1. Recalling that rj = 2sj — x, 
ri2 = 2(si — S2), Pj = 2(x — Sj),j = 1,2, we employ the following expansions of the 
exponentials in Eq. (|22|) 



An 



(23) 



Air 



Im 



(24) 



-a|2s-x|— 6|x-s| 



S — X 



Im 



(25) 



-a|si-S2| 



EA^''^(^i,52)ii:;,(J0>1n^(S2), 



Si — S2 S1S2 



Im 



(26) 



~a|si--S2| 



Air 



Y.BT\s,,S2)Y:^{h)YUs2), 

■5l"52 Ira 



(27) 



where the yj^'s are the usual spherical harmonics. Using Eqs. (^) — (|27D in Eq. ( [2^ ) we 

get 

L 

E 

«=0 



^ Jo Jo pSo Jo 



X 



h{x)G\'''\s„x)Gr>{s2,x)Br'"{s,,S2) - Jr'isux) 



yic,d) 



PS2 

(27n), 



j(a,b). 



X GS^''^)(.2,a:)i?r"^(^i,^2) + -Gr'{s,,x)Gr'{s2,x) 



>(27n)/ 



^(a,fe) , 



-f(c,rf)/ 



X A^^"\s„S2) - -Gr'{si,x)Kr'{s2,x)Br>{suS2) 



,(a,6). 



>(27n)/ 



(28) 



where the /-sum is truncated at Z = L. This procedure avoids complicated angular integra- 
tions involving Si, S2 and x. The matrix element takes a simple form requiring straightfor- 
ward numerical computation of certain radial integrals only. The functions Gi, Ji, Ki etc. 
are easily calculated using Eqs. ( [23| ) — (|27|): 



G\"''''\s,x) = — I duPi{u)e 
2 j-i 



sx 



'a|2s— x|— fe[x— s| 



(29) 



where Pi{u) is the usual Legendre polynomial and u is the cosine of the angle between s and 
X. The integrals ([T8| ) and ( |19|) can be evaluated similarly. For example 



fOO poo roo ^ 

Jo Jo Jo 

X h{x)G\''^\sux)G\''''\s2,x)Bi'''-\suS2) - Jt'^\sux) 



>(2'7'mn ) / 



X A?^-\s„S2) - -Gr\s„x)Kr\s„x)Br->{s„S2] 



)(27mn) / 



(30) 



where e = X + 6m + fJ^n, f = + + I, g = X + 5n + l^m, h = 2am + 27^ + 1 and 



III. NUMERICAL RESULT 

We tested the convergence of the integrals by varying the number of integration points 
in the x, Si and S2 integrals in Eqs. (^) and (|30|) and the u integral in Eq. (|29|) . The 
X integration was relatively easy and 20 Gauss-Legendre quadrature points appropriately 
distributed between and 16 were enough for convergence. In the evaluation of integrals 
of type (p9D 40 Gauss-Legendre quadrature points were sufficient for adequate convergence. 
The convergence in the numerical integration over si and S2 was achieved with 300 Gauss- 
Legendre quadrature points between and 12. The maximum value of / in the sum in Eqs. 
(p8|) and (|30D, L, is taken to be 7 which is sufficient for obtaining the convergence with the 
partial- wave expansions ( p3D — (pT]). 

We find that a judicial choice of the parameters in Eq. (|T3p is needed for convergence. 
The present method does not provide a bound on the result. Consequently, the method 
could lead to a wrong scattering length if an inappropriate (incomplete) basis set is chosen. 
After some experimentation we find that for good convergence the nonlinear parameters 
6n and a„ should be taken to have both positive and negative values and f3n should have 
progressively increasing values till about 1.5. If no care is taken in choosing the parame- 
ters a large number of functions could be necessary for obtaining convergence. The results 



reported in this work are obtained with the following parameters for the functions fn,n ~ 
1,...,14: {Sn,an,(3n,jn,l^n} = {-0.5,-0.25,0.3,0.01,0.02}, {-0.7,-0.25,0.5,0.04,0.02}, 
{-0.7, -0.25, 0.7, 0.03, 0.06}, {-0.4, -0.1, 0.6, 0.2, 0.2}, {-0.2, 0.1, 0.8, 0.2, 0.2}, 
{0.4, -0.2, 0.6, 0.3, 0.3}, {-0.2, -0.1, 0.7, 0.4, 0.3}, {0.3, 0.2, 0.8, 0.3, 0.4}, 

{0.2, 0.2, 1, 0.4, 0.4}, {0.3, 0.2, 1.2, 0.5, 0.5}, {-0.5, 0.2, 1.3, 0.6, 0.6}, {-0.2, 0.1, 1.4, 0.7, 0.7}, 
{0.3,0,1.5,0.8,0.8}, {-0.2,-0.1,1.6,0.9,0.9}. By employing a suitably chosen set of the 
parameters we have kept the number of functions to a minimum. 

Table I: Ps-He scattering length in au for different L and N. 



N 


L = 


L = 2 


L = 3 


L 


= 4 


L 


= 5 


L 


= 6 


L = 7 


1 


2.056 


1.867 


1.782 


1 


721 


1 


681 


1 


655 


1.638 


3 


2.418 


-1.662 


28.016 


4 


151 


3 


061 


2 


839 


2.933 


5 


-13.653 


1.563 


1.234 


1 


128 


1 


135 


1 


168 


1.190 


6 


0.712 


0.982 


0.782 





637 





650 





769 


0.878 


7 


3.372 


0.983 


0.792 





694 





727 





824 


0.910 


8 


-1.657 


1.124 


0.944 





877 





907 





971 


1.023 


9 


6.283 


0.976 


0.832 





801 





841 





897 


0.943 


10 


1.332 


1.123 


0.941 





897 





929 





981 


1.023 


11 


1.225 


1.112 


0.945 





886 





909 





963 


1.011 


12 


0.756 


1.468 


0.995 





918 





931 





964 


0.977 


13 


1.229 


1.197 


1.008 





944 





970 


1 


028 


1.022 


14 


1.061 


1.249 


1.060 





980 





983 


1 


026 


1.019 



In Table I we show the convergence pattern of the present calculation w.r.t. the number of 
partial waves L and basis functions used in the calculation. The convergence is satisfactory 
considering that we are dealing with a complicated five-body problem. However, as the 
present calculation does not produce a bound on the result, the convergence is not monotonic 
with increasing N. The final result of the present calculation is that for = 14 and L = 7: 
a — 1.02 au. Although it is difficult to provide a quantitative measure of convergence, from 



the fluctuation of tliis result for large and L we believe the error in our result to be less 
than 10%, so that the final Ps-He scattering length is taken as a = (1.0 ± 0.1) au. The 
results for large and L reported in Table I all lie in this domain. 

The maximum number of functions (A^ = 14) used in this calculation is also pretty 
small, compared to those used in different Kohn-type variational calculations for electron- 
hydrogen (A^ = 56) positron- hydrogen (A^ < 286) and positron-hehum (A^ < 502) 
||32| scattering. Because of the explicit appearance of the Green's function, the present 
basis-set approach is similar to the Schwinger variational method. Using the Schwinger 
method, convergent results for electron-hydrogen and positron-hydrogen scattering 



have been obtained with a relatively small basis set (A^ ~ 10). These suggest a more rapid 
convergence in these problems with a Schwinger-type method. 

IV. SUMMARY AND DISCUSSION 

To summarize, we have formulated a basis-set calculational scheme for S-wave Ps-He 
elastic scattering below the lowest inelastic threshold using a variational expression for the 
transition matrix. We illustrate the method numerically by calculating the scattering length 
in the electronic doublet state: a = (1.0 ± 0.1) au. This corresponds to a zero-energy 
cross section of (4.0 ± 0.8)7raQ in reasonable agreement with a model calculation by this 
author (3.347raQ) and the experiment of Skalsey et al. [(2.6 ± 0.5)7raQ at 0.9 eV] 0. 
This calculation as well as our previous studies of Ps-He scattering using a model exchange 
potential [0,^] possibly consolidate the experimental result of Skalsey et al. However, these 



low-energy Ps-He elastic scattering cross sections are in disagreement with other experiments 
by Nagashima et al. [(13 ± 4)7ra^ at 0.15 eV] [g, by Canter et al. (8.477rag) [||, and by 
Coleman et al. (Qytoq) [P|, as well as with conventional static-exchange calculations of Refs. 
; TT| , [T2|JT5| , p!7[| (~ 147raQ) and a model potential calculation of Ref. |[T^ (7.737raQ). 



As the effective interaction for elastic scattering between Ps and He is repulsive in nature. 



a smaller scattering length as obtained in this study and in Refs. [0,^ would imply a weaker 



effective Ps-He interaction. This would allow the Ps atom to come closer to He and would 



lead ||25| to a large pickoff quenching rate and a large ^^gff (~ 0.11) in agreement with 
experiment P,|l^. The conventional close-coupling [Q, R-matrix [§] and static-exchange 
|lTTl , p!5| - |T7|] models yielded a much too large scattering length corresponding to a stronger 
repulsion between Ps and He. Consequently, these models led to a much too small ^Z^q 
(~ 0.04) |[T5|- |T7|^ in disagreement with experiment P,[l^. This is addressed in detail in 
Ref. ||2^ where we established a correlation between the different scattering lengths and the 
corresponding ^Z^q. This correlation suggests that a small Ps-He scattering length as in 
this work is consistent with the large experimental ^Z^q. 

Although we have used a simple wave function for He in this complex five-body calcula- 
tion we do not believe that the use of a more refined He wave function would substantially 
change our findings and conclusions. However, independent calculations and accurate ex- 
periments at low energies are welcome for a satisfactory resolution of this controversy. 

The work is supported in part by the Conselho Nacional de Desenvolvimento - Cientifico 
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Figure Caption: 

Figure number 1. Different position vectors for tlie Ps-He system w.r.t. tlie massive 
alplia particle at tfie origin in arrangement 1 witfi electrons 2 and 3 forming He and 1 
forming Ps. The arrows on the electrons indicate the orientations of spin — up and down. 
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